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The spin dynamics of rare earth ions (Er, Nd, Sm and Gd) 
in heavily doped phosphate glasses have been investigated us- 
ing 31 P NMR. Correlation times in the range 10 -2 to 10~ 12 s 
have been obtained for temperatures between 4 K and 100 K. 
No evidence of spin-spin coupling between the ions has been 
found and spin relaxation occurs via conventional phonon pro- 
cesses, including the Orbach process. 

A model involving inhomogeneous broadening of the NMR 
resonance lines, with distinct sample regions corresponding 
to the presence or absence of nuclear spin diffusion, has been 
used in extracting the electron spin correlation times from the 
NMR measurements. 

PACS: 75.50.Kj. 76.30.Kg. 76.60.-k. 



I. INTRODUCTION 

Novel magnetic and magneto-optic phenomena, of in- 
terest for optoelectronic applications, have been found in 
the rare earth (R) metaphosphate glasses (REMG) with 
compositions in the vicinity (R203)o.2s(P205)o.75- For ex- 
ample, these paramagnetic materials exhibit the largest 
known magnetic contributions to the low temperature 
specific heats in oxide glasses 1 . X-ray diffraction and EX- 
AFS studies 2-4 have shown that the structure of REMG 
comprises a 3-D network of corner-linked PO4 tetrahe- 
dra, with the rare-earth ions, which in several instances 
are known to be trivalent R 3+ , occupying sites within the 
PO4 skeleton. 

In recent work 5 , we have shown that it is possible to 
study the spin dynamics of rare earth ions in metaphos- 
phate glasses using 31 P NMR as a probe. The present 
investigation has extended the NMR measurements to a 
number of REMG systems containing Er, Nd, Sm and Gd 
ions. For some of these systems, La or Y ions have been 
used as a buffer in order to lower the concentration of 
magnetic ions while preserving the metaphosphate glass 
structure. Information on the dynamics of the rare earth 
ions has been obtained using the model developed previ- 
ously 5 . 



II. THEORY 

The local fields set up by paramagnetic ions lead to 
nuclear resonance frequency shifts and the establishment 
of an exclusion barrier of radius p Cl inside which spins 
are excluded from the observed NMR signal, together 
with a diffusion barrier of radius bo, inside which nu- 
clear spin diffusion does not operate 5 . In highly param- 
agnetic systems of the type studied here, nuclear spin sys- 
tems are thus composed of three distinct spatial regions 
corresponding to spin diffusive, non-diffusive and exclu- 
sion zones. The theory of nuclear spin lattice relaxation 
(NSLR) due to localized, isolated paramagnetic centres 
in magnetically dilute crystals 6-11 has been extended to 
account for magnetically concentrated systems 5 . In such 
systems, nuclear spin diffusion may be inoperative in ap- 
preciable fractions of sample volume, and nuclei may re- 
lax to multiple paramagnetic sites. The average relax- 
ation rate for nuclei in diffusive regions is given by 
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D is the nuclear spin diffusion coefficient, n s is the mag- 
netic ion concentration and tu n is the nuclear Larmor 
frequency 5 . 

The net electronic relaxation rate for ions with a crys- 
tal field splitting includes direct, Raman and Orbach re- 
laxation rates 12-14 : 
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C d , Cr and Co are relaxation constants for the direct 
, Raman and Orbach processes, respectively. A is the 
crystal field splitting. 
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In Gd 3+ doped glasses, there is no crystal field split- 
ting, and the ground state is a multiplet, so the net re- 
laxation rate may be expressed as 
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Equations f| and || assume a Debye phonon density 
of states. This is consistent with recent work 15 , includ- 
ing the soft potential model approach, which suggests 
that vibrational modes in these systems are phonon-likc. 
Fractons 17 ' 18 and 'excess' modes 15 involving low-energy 
excitations do not appear to play a role in the relaxation 
processes in these systems. 



III. EXPERIMENTAL DETAILS 

The samples were 

prepared from the melt at Bath University with origi- 
nal stoichiometry 24 (R 2 3 ) 2: (X2 03).25-x(P2 5 ).7 5 , and 
x in the range 0.01 to 0.25 (nominally), for R = Er 3+ , 
Nd 3+ , Sm 3+ and Gd 3+ , and X = Y 3 + or La 3+ . 

31 P NMR measurements at 8.5MHz (0.49 T) and 
19.25MHz (1.1 T) were carried out using a pulsed 
NMR spectrometer and a Varian water-cooled electro- 
magnet. Preliminary EPR measurements were carried 
out on Er and Gd REMG using a commercial X-band 
Bruker ESP380E pulsed EPR spectrometer and an Ox- 
ford continuous-flow helium cryostat to achieve low tem- 
peratures. 

Nuclear resonance lineshapes have been characterised 
in terms of the moments of the resonance line deter- 
mined from the Fourier Transform frequency-domain 
NMR spin-echo waveforms for samples with large inho- 
mogeneous broadening. 

The average nuclear relaxation rates for diffusive re- 
gions were extracted from the magnetization recovery 
data as a function of time r using 5 
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where v is a scaling factor, introduced because the de- 
gree of nuclear saturation is not known precisely, and / 
is the diffusive fraction of nuclei. The first and second 
terms describe NSLR in diffusive and non-diffusive re- 
gions, respectively. The relaxation rate for non-diffusive 
regions has been expressed in the form T[ = \T\ for con- 
venience, where A is the ratio of the average relaxation 
rates for nuclei in diffusive and non-diffusive regions. 



IV. RESULTS AND DISCUSSION 

We were unable to observe EPR free induction decay 
or spin-echo signals in the samples investigated, prob- 
ably due to extremely short spin-spin relaxation times 



in these systems. Direct measurements of the electron 
spin-lattice relaxation time using pulsed EPR techniques 
have, therefore, not been possible. 

Nuclear resonance lineshapes have been characterised 
in terms of the moments of the resonance line, using 
Auj =< M% >2 , and the moment ratio, R=< M4 > / < 
A/2 > 2 . The NMR resonance line is inhomogeneously 
broadened due to coupling of the nuclear and paramag- 
netic ion moments. Fig. la shows plots of the inhomo- 
geneouosly broadened 31 P nuclear resonance linewidth 
and moment ratio R versus inverse temperature for a 1% 
Er REMG at B = 1.1 T. The linewidth increases with 
decreasing temperature, eventually reaching a plateau at 
approximately 20 K. The ratio R is close to 3 at low tem- 
peratures, suggesting a characteristically Gaussian line- 
shape. The natural dipolar linewidth is estimated to be 
approximately 3.8 kHz, so there is appreciable inhomo- 
geneous broadening of the linewidth even at the highest 
temperatures at which measurements were made. Fol- 
lowing the onset of motional narrowing, R increases to 
around 5, implying a change in lineshape. This behaviour 
is typical of all the glasses studied. Fig. lb shows plots 
of linewidth versus inverse temperature for various Er 
REMG samples . For each sample, the linewidth in- 
creases while the signal amplitude decreases with decreas- 
ing temperature, eventually becoming unobservably weak 
over a small temperature range. The signal-to-noise ratio 
increases with increasing resonance frequency so that, for 
a given sample, the temperature range over which the line 
can be observed is somewhat larger for a larger applied 
field. 

It appears that the maximum observed linewidth de- 
pends on the magnetic ion concentration, the particular 
ion present and the magnetic field. In all cases, the NMR 
signal became immeasurably broad over a very small tem- 
perature range at sufficiently low temperatures. The rel- 
atively low maximum observed linewidth (Aw < 25 kHz 
) and the difference in maxima for different experiments 
strongly suggest that the abrupt loss of signal at low 
temperatures is not associated with limited spectrome- 
ter bandwidth. This loss of signal appears to be linked 
to the overlapping of diffusion barriers when two or more 
ions contribute to the line-broadening process for a ma- 
jority of nuclear spins in all regions of the sample. 

The 31 P spin-spin relaxation time, T2, determined from 
the spin-echo envelope, is 260 /is. For the Er 3+ , Nd 3+ 
and Sm 3+ ion doped glasses, the electron correlation time 
approaches the nuclear T2 value at temperatures in the 
vicinity of 5 K. For Gd 3+ , the temperature at which r e 
becomes comparable to T2 is approximately 50 K. When 
T e ^ T2, we expect the diffusion and exclusion barrier 
radii to reach their maximum values. 

Fig. 2a shows plots of 31 P relaxation rates T x ~ ver- 
sus T- 1 for the 1% Nd, 6% Nd and metaphosphate Nd 
REMG for B = 4.9 kG. The peak in the relaxation rate 
for the 1% Nd data corresponds to uj n T e — 1. The tem- 
perature range over which measurements could be made 
on the 6% Nd and metaphosphate Nd REMG before the 
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line broadened dramatically rendering the signal unob- 
servable, was insufficient to observe the peak. It is easily 
seen that the nuclear relaxation rate is proportional to 
the magnetic ion concentration, but it appears that the 
metaphosphate Nd REMG contains a magnetic ion con- 
centration lower than the nominal 25%. The data have 
been analysed using Eq.g [DL] and Eq| [RD]. The curves 
represent the best fits to the data of the [RD] and [DL] 
expressions. The Nd 3+ data are well-described by the 
[RD] expression above 20 K, and follow the [DL] expres- 
sion at lower temperatures. In Nd 3+ , the contributions 
to electronic relaxation from the Raman and direct pro- 
cesses are negligible, and the data were analysed taking 
only the Orbach process into account for the temperature 
range covered in these experiments. All three Nd data 
sets were fitted using the same Orbach relaxation param- 
eters. Both the lineshape data and the relaxation data 
provide no evidence that spin-spin interactions between 
the rare earth ions are important over the temperature 
range covered in the present experiment. 

Fig. 2b shows plots of 31 P relaxation rates ver- 
sus T" 1 for the 1% Nd, 1% Er, and 20% Sm REMG. 
The nuclear relaxation rates for the Nd and Er REMG 
are similar at high temperatures. The peak in the Nd 
data occurs at a lower temperature than in the Er and 
Sm data. This suggests faster electronic relaxation in 
the Nd 3+ ions, which results from a smaller crystal field 
splitting than in the other REMG. Nuclear relaxation in 
the 20% Sm REMG is much faster than in the Er and Nd 
REMG because of the much higher dopant concentration. 

For the 1% Gd REMG, electronic relaxation is much 
less efficient than in the other REMG due to the absence 
of the Orbach process. The condition uj n r e — 1 is not sat- 
isfied within the temperature range, and all the data are 
on the low temperature side of the relaxation peak (Fig. 
2c). In Gd 3+ , the whole data set is well described by the 
[DL] expression. The direct relaxation process dominates 
relaxation below 20 K, while the Raman process is most 
important at higher temperatures. Electronic relaxation 
parameters for Nd 3+ , Er 3+ , and Sm 3+ ions appear in Ta- 
ble 1. For Gd 3+ ions, the Orbach process is not operative 
and direct and Raman processes determine the ionic re- 
laxation time. We obtain Cd = 8.1(2) x 10~ 2 (s.K) and 
C R = 2.9(2) xlO 5 (s.K 5 ) for the Gd 3+ REMG. 

The parameters listed in Table 1, together with Cr 
and Cd for Gd 3+ , have been used to generate plots of the 
electronic relaxation rate r e _1 vs. T _1 for the magnetic 
ions studied (Fig. 3). It can be seen that the relaxation 
rates for ions with a crystal field splitting are much higher 
than for the Gd 3+ ions. As noted above, it is clear from 
Fig. 2b that the relaxation behaviour of the 1% Nd and 
1% Er glasses is very similar at higher temperatures. Fig. 
3 shows that a crossover in electron correlation times for 
these two systems occurs at around 15 K. 

The quality of the fits of the relaxation expressions to 
the experimental data and the consistency of the fitting 
parameters between data sets provide strong support for 
the model as well as for the use of electronic relaxation 



expressions assuming a phonon density of states. 



V. CONCLUSION 

Measurements of the 31 P spin- lattice relaxation time 
in metaphosphate glasses doped with rare earth ions have 
provided information on the spin dynamics of these ions 
over a wide range of temperature. For Er, Nd and Sm 
ions, the observed behaviour with temperature is simi- 
lar, with the spin correlation time being determined by 
a two-phonon Orbach process. The local environment of 
the ions is a distorted octahedron of oxygen atoms which 
give rise to crystal field splittings. Estimates of the split- 
ting are provided by the present work for the ions men- 
tioned above. In the case of Gd ions, where no crystal 
field splitting occurs, the paramagnetic ion correlation 
time is determined by direct and Raman processes. The 
approach we have established for determining details of 
paramagnetic ion behaviour in these systems should be 
widely applicable to similar systems. 
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Co (s.K 5 ) 


A(K) 


Nd 3+ 


2.7(3) xl0~ Y 


88(5) 


Er 3 + 


1.6(3) xl0~ 7 


102(5) 


Sm 3 + 


1.1(3) xl0~ 7 


103(5) 



TABLE I. Electronic relaxation parameters for magnetic 
rare earth ions obtained from 31 P relaxation data. The quan- 
tities are defined in Eqs. (4) and (5) 



FIG. 1. (a)Inhomogeneously broadened nuclear resonance 
linewidth (dots) and moment ratios (squares) for a 1% Er 
REMG with B = 1.1 T. (b)Inhomogeneously broadened nu- 
clear resonance linewidths for 1% Er and 5% Er REMG with 
B = 0.49 T and B = 1.1 T. 



FIG. 2. (a)Plots of T^ 1 versus T" 1 for 1% Nd, 6% Nd 
and Metaphosphate Nd REMG. (b)Plots of T" 1 versus T" 1 
for 1% Nd, 1% Er and 20% Sm. (c) A plot of T^ 1 versus 
T _1 for 1% Gd REMG. The curves are the best fits to the 
experimental data of the [RD] and [DL] expressions. 

FIG. 3. Plots of the electronic relaxation rate extracted 
from nuclear relaxation data vs. inverse temperature for the 
magnetic rare earth ions studied. 
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